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The  presence  of  long-range  quantum  spin  correlations  un¬ 
derlies  a  variety  of  physical  phenomena  in  condensed  mat¬ 
ter  systems,  potentially  including  high-temperature  supercon¬ 
ductivity  [1,  2].  However,  many  properties  of  exotic  strongly 
correlated  spin  systems  (e.g.,  spin  liquids)  have  proved  diffi¬ 
cult  to  study,  in  part  because  calculations  involving  /V-body 
entanglement  become  intractable  for  as  few  as  N  ~  30  par¬ 
ticles  [3].  Feynman  divined  that  a  quantum  simulator  —  a 
special-purpose  "analog"  processor  built  using  quantum  par¬ 
ticles  (qubits)  —  would  be  inherently  adept  at  such  problems 
[4,  5].  In  the  context  of  quantum  magnetism,  a  number  of  ex¬ 
periments  have  demonstrated  the  feasibility  of  this  approach 
[6-14].  However,  simulations  of  quantum  magnetism  allow¬ 
ing  controlled,  tunable  interactions  between  spins  localized 
on  2D  and  3D  lattices  of  more  than  a  few  10’s  of  qubits  have 
yet  to  be  demonstrated,  owing  in  part  to  the  technical  chal¬ 
lenge  of  realizing  large-scale  qubit  arrays.  Here  we  demon¬ 
strate  a  variable -range  Ising-type  spin-spin  interaction  Jjj  on 
a  naturally  occurring  2D  triangular  crystal  lattice  of  hundreds 
of  spin- 1/2  particles  (9Be+  ions  stored  in  a  Penning  trap),  a 
computationally  relevant  scale  more  than  an  order  of  magni¬ 
tude  larger  than  existing  experiments.  We  show  that  a  spin- 
dependent  optical  dipole  force  can  produce  an  antiferromag¬ 
netic  interaction ./(  /  °c  where  a  is  tunable  over  0  <  a  <  3; 
dj j  is  the  distance  between  spin  pairs.  These  power-laws  cor¬ 
respond  physically  to  infinite-range  (a  =  0),  Coulomb-like 
(a  =  1),  monopole -dipole  (a  =  2)  and  dipole-dipole  ( a  =  3) 
couplings.  Experimentally,  we  demonstrate  excellent  agree¬ 
ment  with  theory  for  0.05  <  a  <  1.4.  This  demonstration  cou¬ 
pled  with  the  high  spin-count,  excellent  quantum  control  and 
low  technical  complexity  of  the  Penning  trap  brings  within 
reach  simulation  of  interesting  and  otherwise  computationally 
intractable  problems  in  quantum  magnetism. 


A  challenge  in  condensed  matter  physics  is  the  fact  that 
many  quantum-magnetic  interactions  cannot  currently  be 
modeled  in  a  meaningful  way.  A  canonical  example  is  the 
spin-liquid  postulated  by  Anderson  [1].  He  suggested  this  ex¬ 
otic  state  would  arise  in  a  collection  of  spin- 1/2  particles  resid¬ 
ing  on  a  triangular  lattice  and  coupled  by  a  nearest-neighbor 
antiferromagnetic  Heisenberg  interaction.  The  spin-liquid’s 
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Figure  1.  The  Penning  trap  confines  hundreds  of  spin-1/2  particles 
(qubits)  on  a  two-dimensional  (2D)  triangular  lattice.  Each  qubit  is 
the  valence  electron  spin  of  a  yBe+  ion.  (lower)  A  Penning  trap  con¬ 
fines  ions  by  use  of  a  combination  of  static  electric  and  magnetic 
fields.  The  trap  parameters  are  configured  so  that  laser-cooled  ions 
form  a  triangular  2D  crystal.  A  general  spin-spin  interaction  Hi  is 
generated  by  a  spin-dependent  excitation  of  the  transverse  (along  z) 
motional  modes  of  the  ion  crystal.  This  coupling  is  implemented 
with  an  optical  dipole  force  (ODF)  due  to  a  pair  of  off-resonance 
laser  beams  (left  side)  with  angular  separation  Qr  and  difference  fre¬ 
quency  gR.  Microwaves  at  124  GHz  are  directed  to  the  ions  with  a 
horn  and  permit  global  spin  rotations  He-  (upper)  A  representative 
top-view  resonance-fluorescence  image  showing  the  center  region  of 
an  ion  crystal  captured  in  the  ions’  rest-frame;  in  the  lab-frame  the 
ions  rotate  a  at  co;-=  lit  x  43.8  kHz  [21],  Fluorescence  is  an  indica¬ 
tion  of  the  qubit  spin-state  (Jf)  bright,  |4_)  dark);  here,  the  ions  are  in 
the  state  ||).  The  lattice  constant  is  do  ~  20  pm. 


ground  state  is  massively  degenerate,  owing  to  spin  frustra¬ 
tion,  and  is  expected  to  exhibit  unusual  behaviors  including 
phase  transitions  at  zero  temperature  driven  by  quantum  fluc¬ 
tuations  [15].  However,  despite  recent  advances  [16,  17]  a 
detailed  understanding  of  large-scale  frustration  in  solids  re¬ 
mains  elusive  [2,  18-20]. 
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Atomic  physicists  have  recently  entered  the  fray,  providing 
a  bottom-up  approach  by  engineering  the  relevant  spin  inter¬ 
actions  in  quantum  simulators  [5,  22,  23].  The  necessary  ex¬ 
perimental  capabilities  were  first  demonstrated  in  the  context 
of  atomic  clocks:  laser  cooling,  deterministic  spin  localiza¬ 
tion,  precise  spin-state  quantum  control,  high-fidelity  readout, 
and  engineered  spin-spin  coupling  (e.g.,  [24]).  In  the  domain 
of  quantum  magnetism,  this  tool  set  permits  control  of  param¬ 
eters  commonly  viewed  as  immutable  in  natural  solid  systems, 
e.g.,  lattice  spacing  and  geometry,  and  spin-spin  interaction 
strength  and  range. 

Initial  simulations  of  quantum  Ising  and  Heisenberg  inter¬ 
actions  with  localized  spins  were  done  with  neutral  atoms  in 
optical  lattices  [6, 11],  atomic  ions  in  Paul  traps  [9, 10, 13, 14], 
and  photons  [12],  This  work  involved  a  domain  readily  cal¬ 
culable  on  a  classical  computer:  interactions  between  N  ~  10 
qubits  localized  in  ID  chains.  The  move  to  quantum  magnetic 
interactions  on  2D  lattices  and  larger,  computationally  rele¬ 
vant  particle  numbers  is  the  crucial  next  step  but  can  require 
significant  technological  development  [25]. 

In  our  Penning  trap  apparatus,  laser-cooled  9Be+  ions  natu¬ 
rally  form  a  stable  2D  Coulomb  crystal  on  a  triangular  lattice 
with  ~  300  spins  (Fig.  1).  Each  ion  is  spin-1/2  system  (qubit) 
over  which  we  exert  high  fidelity  quantum  control  [26].  In 
this  paper  we  demonstrate  the  use  of  a  spin-dependent  optical 
dipole  force  (ODF)  to  engineer  a  continuously-tunable  Ising- 
type  spin-spin  coupling  Jt  j  dl  / .  This  capability,  in  tandem 
with  a  modified  measurement  routine  (e.g.,  by  more  sophisti¬ 
cated  processing  of  images  like  Fig.  1),  is  a  key  advance  to¬ 
ward  useful  simulations  of  quantum  magnetism. 

A  Penning  trap  confines  ions  in  a  static  quadrupolar  electric 
potential  (see  Methods)  and  a  strong,  homogeneous  magnetic 
field  Bqz  (here,  Bq  =  4.46  Tesla).  Axial  trapping  (along  z)  is 
due  to  the  electric  field.  Ion  rotation  at  frequency  cnr  (about 
z)  produces  a  radial  restoring  potential  due  to  the  velocity- 
dependent  Forentz  force  (qv  x  B).  Tuning  the  ratio  of  the  ax¬ 
ial  to  radial  confinement  allows  controlled  formation  of  a  pla¬ 
nar  geometry  and,  after  Doppler  laser  cooling,  the  formation 
of  a  2D  Coulomb  crystal  on  a  triangular  lattice  (see  [21]  and 
Methods).  We  routinely  generate  crystals  with  100  <  N  <  300 
ions,  where  the  valence-electron  spin  state  of  each  ion  serves 
as  a  qubit  [26].  Following  techniques  developed  in  linear  (ID) 
Paul  traps  [27],  spins  confined  in  the  same  trapping  potential 
are  coupled  through  their  shared  motional  degrees  of  freedom. 

Using  well  controlled  external  fields,  we  engineer  spin  in¬ 
teractions  of  the  form 

H n  =  D;  B/j  '  /  j  \ 

Hi  = 

where  oj  is  the  z-Pauli  matrix  for  ion  i.  We  label  the  qubit  spin 
states  |t)  =  | ms  =  +1  /2)  and  |j.)  =  | ms  =  —1/2),  where  ms  is 
the  spin’s  projection  along  the  Bqz  quantizing  field,  so  that 
of  |ti)  =  +  |t/)  and  Iti)  =  ^  |4-/)-  Hb  is  an  interaction  (e.g., 
Bu6a)  generated  by  externally  applied  microwaves  at  124  GHz 
that  couples  equally  to  all  spins  and  permits  global  rotations 
(Fig.  1).  The  interaction  Hi  describes  a  general  coupling  Jjj 


a) 


Figure  2.  Spin-spin  interactions  are  mediated  by  the  ion  crystal’s 
transverse  motional  degrees  of  freedom,  (a)  For  a  2D  crystal  with 
N  =  217  ions  and  to,-  =  2tt  x  45.6  kHz,  we  calculate  the  eigenfunc¬ 
tions  bm  and  eigenfrequencies  to,,,  for  the  N  transverse  motional 
modes  (see  Supplementary  Information).  Plotted  here  are  co,„  and 
b,„  for  the  14  highest  frequency  modes.  Relative  mode  amplitude  is 
indicated  by  color.  The  center-of-mass  motion  (COM)  is  the  highest 
in  frequency  (coi  ~  2n  x  795  kHz)  and  has  no  spatial  variation  in  b,„. 
The  lowest-frequency  mode  (CO217  ~  2ji  x  200  kHz)  exhibits  spatial 
variation  in  bm  at  the  lattice  spacing  length  scale  do  ~  20pm.  (b)  By 
use  of  Eq.  3,  Jjj  is  calculated  explicitly  for  IV  =  217  spins  and  plotted 
as  a  function  of  spin-spin  separation  djj.  For  /j«  —  a>i  <  271  x  1  kHz, 
Hqdf  principally  excites  COM  motion  in  which  all  ions  equally  par¬ 
ticipate:  the  spin-spin  interaction  is  spatially  uniform.  As  the  detun¬ 
ing  is  increased,  modes  of  higher  spatial  frequencies  participate  in  the 
interaction  and  Jjj  develops  a  finite  interaction  length.  We  find  the 
scaling  of  Jjj  with  spin  separation  djj  follows  the  power  law  dff. 
For  /jr  —  (ti  1  2nx  500kHz,  all  transverse  modes  participate  and 
the  spin-spin  coupling  power  law  approaches  a  =  3.  The  solid  lines 
are  power-law  fits  to  the  theory  points.  For  comparison  with  other 
experiments,  the  nearest  neighbor  coupling  (do  =  20 gm)  is  denoted 
by  the  dashed  line,  (c)-(e)  The  power-law  nature  of  Jj  j  is  qualita¬ 
tively  illustrated  for  N  =  19  (larger  N  is  illegible).  Nodes  (spins)  are 
joined  by  lines  colored  in  proportion  to  their  coupling  strength  for 
various  a.  (f)  For  context,  the  graph  for  a  ID  nearest-neighbor  Ising 
interaction,  a  well-known  model  in  quantum  field  theory. 
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between  spins  i  and  j,  distance  djj  apart  [28,  29];  7,  ,/  >  o 
(J,  j  <  0  )  is  an  antiferromagnetic  (ferromagnetic)  coupling. 

We  implement  77/  using  a  spatially  uniform  spin-dependent 
optical  dipole  force  (ODF)  generated  by  a  pair  of  off-resonant 
laser  beams  with  difference  frequency  /jr  (see  Fig.  1  and  Sup¬ 
plementary  Information).  The  ODF  couples  each  ion’s  spin 
to  one  or  more  of  the  N  transverse  (along  z)  motional  modes 
of  the  Coulomb  crystal  by  forcing  coherent  displacements  of 
the  ions  that  in  turn  modify  their  Coulomb  potential  energy 
through  the  interaction 

N 

Hodf  =  (2) 

i 

Here  Fz(t)  =  Co  cos /./^f  is  the  ODF,  Zi  = 
Im=i  +a1me,a'"‘)  is  the  axial  posi¬ 

tion  operator  for  ion  i,  and  /?,■  „,  are  elements  of  the  N 
transverse  eigenfunctions  bm  at  frequencies  co,„  normalized  as 
Lm=t  |^,rn|2  =  LoLi  \bi,m\2  =  1  [28,  29].  The  modes  include 
the  center-of-mass  mode  (<Bi)  as  well  as  an  array  of  higher 
spatial-frequency  modes  that  may  be  derived  from  atomistic 
calculations  (see  Fig.  2)  and  confirmed  by  experimental 
measurement  [30]. 

For  small  coherent  displacements,  where  residual  spin- 
motion  entanglement  can  be  neglected  (see  [29]  and  Meth¬ 
ods),  Hodf  is  equivalent  to  77/  in  Eq.  1:  spins  i  and  j  are  cou¬ 
pled  in  proportion  to  their  spin  states  6)  and  a:j  and  mutual 
participation  in  each  motional  mode  m.  The  coupling  coeffi¬ 
cient  is  given  by  [29] 


772  AT  N 

J. .  -  M  y 
J'  J  -  2- 


2  TiM 


b  i  ji/b  j  jn 

1  4  -  < 


(3) 


These  pairwise  interaction  coefficients  7,  ,■  can  be  calculated 
explicitly  by  use  of  Eq.  3  and  classical  calculations  of  ion 
motional  modes.  We  find  that  the  range  of  interaction  can 
be  modified  by  detuning  away  from  the  COM  mode  as  dis¬ 
cussed  in  Fig.  2.  In  the  limit  /jr  —  ©i  271  x  500  kHz,  all 
modes  participate  equally  in  the  interaction  and  7,.  ;  dj  j ,  as 
discussed  in  [28].  At  intermediate  detuning,  we  find  a  power- 
law  scaling  of  the  interaction  range  7,.;  where  a  can 

be  tuned  0  <  a  <  3.  That  is,  by  adjusting  the  single  experi¬ 
mental  parameter  /jr  we  can  mimic  a  continuum  of  physical 
couplings  including  important  special  cases:  a  =  0  is  infinite- 
range,  a  =  1  is  monopole-monopole  (Coulomb-like),  a  =  2  is 
monopole -dipole  and  a  =  3  is  dipole-dipole.  Note  that  a  =  0 
results  in  the  so-called  7r  interaction  that  gives  rise  to  spin¬ 
squeezing  and  is  used  in  quantum  logic  gates  (see  Supple¬ 
mentary  Information)  [27].  In  addition,  tuning /jr  also  permits 
access  to  both  antiferromagnetic  (AFM,  /jr  >  COi)  and  ferro¬ 
magnetic  (FM,  (02  <  COi)  couplings  [13]. 

Experimentally,  we  demonstrate  a  tunable-range  Ising  in¬ 
teraction  by  observing  a  global  spin  precession  under  appli¬ 
cation  of  77/.  We  compare  experimental  data  with  the  mean 
field  (MF)  prediction  that  the  influence  of  77/  on  spin  j  can  be 
modeled  as  an  excess  magnetic  field  in  the  ^-direction  Bj  due 
to  the  remaining  A  —  1  spins,  Bj  =  ^  ;  Jij  )  (see  Sup¬ 

plementary  Information).  For  a  general  qubit  superposition 


Figure  3.  Benchmarking  the  2D  Ising  interaction,  (a)  Spin- 
precession  benchmarking  sequence  for  77/.  The  spins  are  prepared 
at  the  outset  in  If)  (a  ferromagnetic  state).  The  spin-spin  interac¬ 
tion  77/  is  present  when  the  ODF  laser  beams  are  on.  We  choose 
fiR  —  CO/  =  n  x  2n/%arm  so  that  for  small  detunings  from  the  COM 
mode  (COi),  the  spins  are  decoupled  from  the  motion  at  the  end  of 
the  drive,  (b)  Evolution  of  a  single  spin  prior  to  the  application  of  the 
spin-echo  Tt-pulse.  c,d)  Plots  exhibiting  spin  precession  (az)  due  to 
77/  as  a  function  of  initial  tipping  angle  0] .  The  error  bars  are  statisti¬ 
cal  (s.d.,  n  ~  200).  The  plots  show  typical  experimental  data  (black) 
and  single-parameter  fits  to  Eq.4.  For  an  antiferromagnetic  coupling 
(AFM)  and  xarm  =  250/js,  /jr  —  COi  =  271  x  4.0  kHz,  7r  =  1.4W-cm~- 
we  obtain  J/Ir  =  271  x  25  Hz  •  W~2  •  cm4  (yellow  fit).  Longer  drive 
periods  and  higher  laser  intensity  7r  yield  a  larger  precession.  For 
"Wm  =  350/rs,  /ar  —  ct)[  =  271  x  2.9  kHz,  7r  =  1.9W-cm~2  we  ob¬ 
tain  7/7J  =  271  x  55  Hz  •  W~2  ■  cm4  (red  fit).  The  data  in  this  plot 
is  typical  of  the  experiments  conducted  for  benchmarking,  (d)  For  a 
much  stronger  interaction,  Eq.  4  cannot  be  used  to  obtain  7  since  the 
MF  assumption  is  no  longer  valid  (see  Supplementary  Information). 
Also,  here  we  used  a  small  negative  detuning  (CO2  <C  jUr  <  (0i)  which 
gives  a  long-range  ferromagnetic  (FM)  interaction.  For  these  exper¬ 
iments  we  set  (0 ,  =  271  x  45.6  kHz.  (e)  Benchmarking  results  for  an 
ion-crystal  with  N  =  206  ±10  ions.  Each  point  is  generated  by  mea¬ 
suring  7  as  in  c)  and  measuring  the  laser  beam  intensity  at  the  ions 
Ir.  The  vertical  error  bars  are  dominated  by  uncertainty  in  Ir  (see 
Supplementary  Information).  The  solid  line  (red)  is  the  prediction 
of  MF  theory  that  accounts  for  couplings  to  all  N  transverse  modes; 
there  are  no  free  parameters.  The  line’s  breadth  reflects  experimental 
uncertainty  in  the  angle  0r  =  4.8°  ±0.25.  The  MF  prediction  for  the 
average  power-law  scaling  coefficient  a  is  drawn  in  green  (right  axis, 
linear  scale). 
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state,  Bj  gives  rise  to  spin  precession  about  z  in  excess  of  that 
expected  due  to  simple  Larmor  precession.  The  experiment 
sequence  shown  in  Fig.  3a  measures  this  excess  precession, 
averaged  over  all  spins  in  the  crystal.  At  the  outset,  each  spin 
is  prepared  in  jf),  followed  by  a  rotation  about  the  x-axis  by 
an  angle  0i .  Hi  is  applied  during  the  arms  of  a  spin-echo;  pre¬ 
cession  (a?)  coherently  adds  over  both  spin-echo  intervals. 
The  final  7t/2-pulse  maps  precession  out  of  the  initial  vz-plane 
into  excursions  along  z  (above  or  below  the  equatorial  plane 
of  the  Bloch  sphere)  that  are  resolved  by  projective  spin  mea¬ 
surement  along  z. 

As  a  function  of  “tipping-angle”  0i  we  detect  global  state- 
dependent  fluorescence  (|f)  bright,  |j,)  dark)  using  a  photo¬ 
multiplier  tube.  This  measurement  permits  a  systematic  study 
of  MF-induced  spin  precession  averaged  over  all  particles, 

=  2  cos0,  =  2/ cos 0i.  The  prob¬ 

ability  of  detecting  |f)  at  the  end  of  the  sequence  is 

p(lt))  =  l  (1  +exp(-r-2xarm)sin(0i)sin(2/cos(0i)  ■2'Carm)) , 

(4) 

and  a  single-parameter  fit  to  experimental  data  yields  J.  De¬ 
coherence  due  to  spontaneous  emission  is  accounted  for  by 
r  and  is  fixed  by  independent  measurement  of  the  ODF  laser 
beam  intensities  Ir  (r  Ir,  see  Supplementary  Information). 

In  Fig.  3c,d  we  show  representative  measurements  of  ex¬ 
cess  precession  due  to  ///  for  different  values  of  spin  cou¬ 
pling  strength  (determined  by  /j  //,)  and  interaction  dura¬ 
tion  2%arm.  The  excess  spin  precession  varies  periodically  with 
01  (period  n)  and  larger  interaction  strength  yields  greater  pre¬ 
cession,  manifested  in  our  experiment  as  a  larger  amplitude 
modulation  of  P(|f)).  Data  agree  with  Eq.  4  and  allow  direct 
extraction  of  /  for  given  experimental  conditions.  In  Fig.  3e 
we  plot  /,  normalized  by  /J  (Ir  independently  measured),  as 
function  of  detuning  /jr  —  00i  (N  =  206  ±  10  ions).  Using  no 
free  parameters,  we  find  excellent  agreement  with  /  obtained 
by  averaging  over  all  J,j,  where  the  J,j  were  calculated  by 
including  couplings  to  all  N  transverse  modes  (Eq.  3). 

The  MF  interpretation  of  our  benchmarking  measurement 
tolerates  only  weak  spin-spin  correlations.  Therefore,  in 
the  benchmarking  regime  we  apply  a  weak  interaction  (/• 
2z/Inn  <§C  VN /4,  see  Supplementary  Information).  In  a  quan¬ 
tum  simulation,  the  same  interaction  is  applied  at  greater 
power  producing  quantum  spin-spin  correlations.  In  the 
present  configuration  of  our  apparatus,  spontaneous  emission 
due  to  the  ODF  laser  beams  (r  in  Eq.  4)  is  the  dominant  source 
of  decoherence.  With  modest  laser  powers  of  4mW/beam  and 
a  detuning  /jr  —  ( Oj  =  2tt  x  2kHz,  we  obtain  /~  2k  x  0.5  kHz 
and  r//~  0.06.  The  expected  spin-squeezing  (Jr)  due  to  this 
interaction  is  5  dB,  limited  by  spontaneous  emission.  The  ratio 
r //can  be  reduced  (/increased)  by  a  factor  of  50  by  increas¬ 
ing  Qr  to  35°,  a  likely  prerequisite  for  access  to  the  shortest- 
range,  dipole-dipole  coupling  regime  (a  — >  3).  At  present,  ge¬ 
ometric  constraints  limit  0^  <  5°;  we  plan  upgrades  to  our 
apparatus  to  permit  Or  =  35°.  We  also  note  that  relaxation  of 
the  constraint  F-  ~  — Fj,  can  also  reduce  T/J. 

In  summary,  this  work  establishes  the  suitability  of  a  Pen¬ 
ning  trap  apparatus  to  pursue  quantum  simulation  in  a  regime 


inaccessible  to  classical  computation.  Our  approach  is  based 
on  naturally  occurring  2D  Coulomb  (Wigner)  crystals  with 
hundreds  of  ion  qubits,  a  novel  experimental  system  that  does 
not  require  demanding  trap-engineering  efforts.  Experimen¬ 
tally,  we  used  an  optical-dipole  force  to  engineer  a  tunable- 
range  spin-spin  interaction  and  benchmarked  its  interaction 
strength.  Excellent  agreement  is  obtained  with  the  predictions 
of  mean  field  theory  and  atomistic  calculations  that  predict 
a  power-law  antiferromegnetic  spin  coupling  JLj  d\  "  for 
0-05  <a<  1.4. 

With  this  work  as  a  foundation,  we  anticipate  a  variety  of 
future  investigations.  For  example,  simultaneous  application 
of  non-commuting  interactions  Hr  and  Hj  is  expected  to  give 
rise  to  quantum  phase  transitions;  Hj  may  be  antiferromag¬ 
netic  (/jr  >  COj)  or  ferromagnetic  (©2  /jr  <  ©1)).  Geometric 
modifications  to  our  apparatus  will  permit  access  to  larger  0^ 
and  antiferromagnetic  dipole -dipole-type  couplings  (a  — >  3). 
Improved  image  processing  software  will  permit  direct  mea¬ 
surement  of  spin-spin  correlation  functions  using  our  existing 
single-spin-resolving  imaging  system  (Fig.  1). 


METHODS 

In  a  frame  rotating  at  frequency  ©,-,  the  trap  potential  is 

qty(r,z)  =  \m&\(z2  +  P>'2),  (5) 

where  q  is  the  ion  charge,  M  is  the  single-ion  mass  and 
P  =  ©,C0|  2(Ylc  —  ©,-)  —  1/2.  The  9Be+  cyclotron  frequency 
is  Yl,  =  Boq/M  =  2n  x  7.6  MHz  and  the  frequency  of  the 
ions’  harmonic  center-of-mass  (COM)  motion  along  z  is  ©1  = 
2n  x  795  kHz.  Ion  rotation  is  precisely  controlled  with  an  ex¬ 
ternal  rotating  quadrupole  potential  [21].  For  100  <  N  <  300 
we  set  ©r  ~  2k  x  45  kHz  so  that  the  radial  confinement  is  suf¬ 
ficiently  weak  that  a  cloud  of  ions  relaxes  into  a  single  2D 
plane  (P  <C  1).  Upon  Doppler  laser  cooling  the  ions’  motional 
degrees  of  freedom  ( Tqom  ~  1  mK)  [30],  the  ions  naturally 
form  a  2D  Coulomb  crystal  on  a  triangular  lattice,  the  geom¬ 
etry  that  minimizes  the  energy  of  their  mutual  Coulomb  po¬ 
tential  energy.  The  crystal  has  N  transverse  eigenmodes  ©,„ 
with  eigenfunctions  bm ;  the  COM  mode  ©1  is  the  highest  fre¬ 
quency  mode  (see  Fig.  2). 

The  spin-dependent  ODF  is  generated  by  a  pair  of  off- 
resonance  laser  beams  with  angular  separation  Or  ~  4.8°  and 
difference  frequency  ,ur  (see  Fig.  1).  The  result  is  a  traveling 
ID  optical  lattice  of  wavelength  A ,r  =  2ji/|a£|  «  3.7  fjm  whose 
wavefronts  propagate  along  z,  traversing  the  ion  crystal  at  fre¬ 
quency  hr /2k.  Alignment  of  A k  is  crucial  for  proper  spin- 
spin  coupling  (see  Supplementary  Information).  The  lattice’s 
polarization  gradient  induces  a  differential  AC  Stark  shift  on 
the  qubit  states  (a  spin-dependent  force).  We  choose  operat¬ 
ing  conditions  that  give  F/  ps  — F/,  where  //  =  Focos(/j/;f)z. 
For  reference,  if  the  single-beam  intensity  at  the  ions  is  Ir  = 
1 W  •  cm'2,  we  obtain  Fq  ~  1 .4  x  10~23  N. 

Small  coherent  displacements  that  produce  negligible  spin- 
motion  entanglement  (as  required  by  Eq.  3)  are  obtained  for 
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detunings  satisfying 


ti  l/u?  iO]tl  |  >  Fq  \fh (2 nm  1 )  /  (2Af  com ) .  (6) 

This  is  a  more  stringent  criterion  than  that  used  by  others 
[13,  29]  as  it  includes  an  additional  \/N  to  account  for  a 
typical  distribution  of  composite  spin  states.  Moreover,  we 
also  include  a  correction  factor  for  finite  temperature  nm  ~ 
k bT 
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SUPPLEMENTARY  INFORMATION 

A  review  of  ion  confinement  in  Penning  traps  and  discus¬ 
sion  of  a  variety  of  equilibrium  states,  including  2D  Coulomb 
crystals  (ionic  Wigner  crystals),  can  be  found  in  [21,  31,  32], 
A  number  of  authors  have  theoretically  analyzed  and  dis¬ 
cussed  the  prospects  of  using  2D  Coulomb  crystals  for  quan¬ 
tum  information  and  computation  [28,  33-37],  The  engi¬ 
neered  Ising  interaction,  which  we  report  here,  builds  on  our 
previous  experimental  work  using  2D  Coulomb  crystals  for 
high-fidelity  quantum  control  [26,  38-41],  Below  we  discuss 
some  details  of  this  new  capability. 


I.  SPIN  INITIALIZATION,  CONTROL,  AND 
MEASUREMENT 

Reference  [26]  gives  a  detailed  discussion  of  our  spin  ini¬ 
tialization,  control,  and  measurement  capabilities  with  planar 
ion  arrays  in  Penning  traps.  Here  we  briefly  summarize  some 
of  that  discussion,  emphasizing  aspects  relevant  for  the  mea¬ 
surements  reported  here.  Figure  4  shows  the  relevant  9Be+ 
energy  levels.  We  use  the  valence  electron  spin  states  paral¬ 
lel  |t)  =  | mj  =  +j)  and  antiparallel  |t)  =  | mj  =  —  j)  to  the 
applied  magnetic  field  of  the  Penning  trap  as  the  two-level 
system  or  qubit.  In  the  4.46  T  magnetic  field  of  the  trap, 
these  levels  are  split  by  approximately  Do  =  2n  x  124  GHz. 
The  9Be+  nucleus  has  spin  1  =  3/2.  However,  we  optically 
pump  the  nuclear  spin  to  the  mj  =  +3/2  level  [42],  where  it 
remains  throughout  the  duration  of  an  experiment.  The  ions 
are  Doppler  laser-cooled  to  a  temperature  ~  1  mK  [30]  by  a 
1  MHz  linewidth,  313  nm  laser  tuned  approximately  10  MHz 
below  the  |f)  — >  \2P?,/2mJ  —  +3/2)  cycling  transition.  Spins 
in  the  |j,)  state  are  efficiently  optically  pumped  to  the  |f)  state 
by  a  laser  tuned  to  the  |j,)  -A  =  +1/2)  transition. 

The  repump  beam  and  the  main  Doppler  laser  cooling  beam 
are  directed  along  the  magnetic  field  (f-axis).  Powers  are  a 
few  milli-Watt  with  laser  beam  waists  of  ~  1  mm.  In  ad¬ 
dition,  a  weak  Doppler  laser  cooling  beam  (~  40  gm  waist) 
directed  perpendicularly  to  the  z-axis  directly  Doppler  cools 
the  perpendicular  degrees  of  freedom.  A  typical  experimental 
cycle  starts  with  10  ms  to  20  ms  of  combined  Doppler  laser 
cooling  and  repumping.  The  repump  laser  remains  on  for  an¬ 
other  3  ms  after  the  Doppler  cooling  laser  is  turned  off.  The 
fidelity  of  the  |f )  state  preparation  is  estimated  to  be  very  high 
(>99.9%)  [26]. 

Low-phase-noise  microwave  radiation  at  124  GHz  is  used 
to  rotate  the  spins  through  the  magnetic  dipole  interaction 

Hg  =  gUBZiB^t)  ■  (gi/2)  ,  where  B^t)  is  the  applied  mi¬ 
crowave  field  (predominantly  perpendicular  to  z),  g  —  2  is  the 
electron  g-factor,  and  /jb  is  the  Bohr  magnetron.  The  fidelity 
of  a  7t-pulse  was  measured  to  be  greater  than  99.9%  in  a  ran¬ 
dom  benchmarking  experiment  [26],  The  microwave  source 
consisted  of  an  agile  15.5  GHz  source  followed  by  an  ampli¬ 
fier  multiplier  chain  with  150  mW  output  power  at  124  GHz. 
The  15.5  GHz  source  is  obtained  by  mixing,  with  a  single- 
side-band  mixer,  the  output  of  a  15.2  GHz  dielectric  resonator 


oscillator  (DRO)  with  the  output  of  a  300  MHz  direct  digi¬ 
tal  synthesizer  (DDS)  that  is  under  field  programmable  gate 
array  (FPGA)  control.  On/off  switching  of  the  124  GHz  mi¬ 
crowaves  is  done  at  15.5  GHz  before  the  amplifier  multiplier 
chain.  The  microwaves  are  transported  to  the  ions  down  the 
bore  of  the  magnet  with  a  rigid  waveguide  and  directed  onto 
the  ions  with  a  horn  located  between  the  ring  and  endcap  elec¬ 
trodes  of  the  trap.  With  this  arrangement  the  microwave  hard¬ 
ware  does  not  block  optical  access  along  the  magnetic  field 
axis,  enabling  imaging  of  the  ion  resonance  fluorescence  scat¬ 
tered  along  the  magnetic  field  (top-view  image  —  see  III). 
We  obtain  7t-pulses  of  70, us  duration  with  the  150  mW  out¬ 
put  power  of  the  amplifier  multiplier  chain  and  the  setup  de¬ 
scribed  here.  The  measured  spin-echo  coherence  duration  ( T2 ) 
is  ~  100  ms. 

At  the  end  of  an  experimental  sequence  we  turn  on  the 
Doppler  cooling  laser  and  make  a  projective  measurement  of 
the  ion  spin  state  through  state-dependent  resonance  fluores¬ 
cence.  With  the  Doppler  cooling  laser  on,  an  ion  in  the  |f) 
state  scatters  ~  107  photons/s  while  an  ion  in  |j,)  is  dark.  For 
the  spin  precession  measurements  reported  here  we  performed 
a  global  spin-state  detection.  Specifically  we  detected,  with 
f/5  light  collection  and  a  photomultiplier  tube,  the  resonance 
fluorescence  from  all  the  ions  in  a  direction  perpendicular  to 
the  magnetic  field  (the  side-view).  For  detection  periods  of 
~  50  ms  the  detection  fidelity  is  high,  typically  limited  by 
quantum  projection  noise.  Here  we  used  short  detection  pe¬ 
riods  of  ~  500/,/s,  from  which  we  would  detect  ~  1  photon 
for  each  bright  state  |f).  Typically,  each  experiment  was  re¬ 
peated  ~  100  times  and  averaged,  resulting  in  a  few  percent 
uncertainty  due  to  shot  noise  in  the  measurement  of  P(t). 

The  spin-precession  signal  used  to  benchmark  spin-spin 
coupling  in  the  manuscript  relied  on  a  global  spin-state  mea¬ 
surement  via  side-view  fluorescence  collected  on  a  photo¬ 
multiplier  tube.  In  the  future,  we  anticipate  that  time -resolved 
top-view  images  such  as  that  shown  in  Fig.  1  will  be  used  to 
obtain  the  spin  state  of  individual  ions.  As  discussed  in  the 
manuscript,  ion  crystal  rotation  at  0+  is  phase-locked  locked 
to  an  external  oscillator.  We  use  an  imaging  photomultiplier 
tube  to  record  (. x,y,t )  for  each  photon.  Rotating-frame  im¬ 
ages  are  generated  computationally  given  (x,y.t)  and  cor,  a 
technique  established  in  2001  [43].  In  linear  Paul  trap  ex¬ 
periments,  determination  of  ion  spin-state  is  possible  with  as 
few  as  10  photons/ion  [44].  At  present  we  await  the  arrival 
of  a  new  (x.y.t)  detector  system  capable  of  a  detection  rate  of 
5  x  106  Hz.  We  anticipate  this  will  enable  high-fidelity  spin- 
state  measurement  of  a  300-ion  crystal  in  ~  10  ms.  From  a 
suite  of  identical  experiments,  the  spin-spin  correlation  func¬ 
tion  can  be  computed.  We  believe  that  our  ability  to  resolve 
single  ions,  even  in  the  presence  of  rotation  at  to r  indicates 
a  path  forward  in  performing  individually  resolved  measure¬ 
ments  of  fluorescence  correlations  between  ions. 

Previous  measurements  have  elucidated  some  of  the  possi¬ 
ble  limiting  mechanisms.  For  instance,  the  stability  of  crys¬ 
tal  orientation  in  the  rotating  frame  was  studied  for  spheri¬ 
cal  crystals  of  ~  10,000  ions.  In  these  experiments  the  ori¬ 
entation  was  observed  to  precess  uniformly  for  durations  of 
>  10  s,  then  suddenly  slip  by  a  large  angle  before  again  re- 
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suming  slow  precession  [43].  This  “stick-slip”  motion  can 
be  followed  and  easily  corrected.  Figure  la  of  the  manuscript 
was  generated  from  600  s  of  integration  after  stick-slip  correc¬ 
tions.  Other  potential  issues  include  ion  loss  and  background 
gas  collisions.  As  the  trap  depth  is  7§>  1  eV,  background  gas 
collisions  do  not  result  ion  loss.  Collisions  with  background 
hydrogen  generate  BeH+  ( 1  per  6  minutes  for  N  ~  300  9Be+ 
ions)  which  collect  at  the  crystal  perimeter  due  to  centrifu¬ 
gal  separation.  These  effects  do  not  appear  to  limit  top-view 
imaging  fidelity,  as  demonstrated  in  Fig.  la  of  the  manuscript. 


II.  OPTICAL-DIPOLE-FORCE  LASER  SETTINGS 

Figure  5  shows  a  simple  sketch  of  the  optical  dipole  force 
(ODF)  laser  beam  set-up.  As  discussed  below,  the  frequency 
as  well  as  the  beam  polarizations  were  chosen  to  null  the 
AC  Stark  shift  from  an  individual  beam  and  to  produce  a  state- 
dependent  force  that  is  equal  in  magnitude  but  opposite  in  sign 
for  the  |t)  and  |t)  qubit  states  (//  =  —F±  ).  This  reduced  the 
system  sensitivity  to  laser  intensity  fluctuations.  For  exam¬ 
ple,  if  Ef  7^  — C ,  then  the  interaction  induced  by  the  optical 
dipole  force  will  include  terms  linear  in  the  oj’s.  These  terms 
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Figure  4.  Relevant  energy  levels  of  yBe+  at  f?o  =  4.46  T  (not  drawn 
to  scale).  We  show  only  mi  =  +  |  levels  that  are  prepared  experimen¬ 
tally  through  optical  pumping.  The  2S\j2  —  2Pi/2  transition  wave¬ 
length  is  ~  313  nm.  A  resonant  laser  beam  provides  Doppler  laser 
cooling  and  state  discrimination;  a  second  laser  beam  repumps  |4) 
to  |t).  The  optical  dipole  force  (ODF)  interaction  is  due  to  a  pair 
of  beams  (derived  from  the  same  laser)  with  relative  detuning 
The  qubit  splitting  AEg/h  ~  271  x  124  GHz.  A  low-phase-noise  mi¬ 
crowave  source  at  124  GHz  provides  full  global  control  over  spins. 


Figure  5.  Sketch  of  optical  dipole  force  (ODF)  laser  beam  setup, 
(a)  The  ODF  laser  beams  lie  in  the  y-z  plane  at  angles  ±9^/2  with 
respect  to  the  y-axis.  (b)  View  looking  in  the  —  y  direction.  The 
beams  are  linearly  polarized  but  with  different  polarization  angles 
relative  to  vertical  polarization. 

can  be  canceled  with  spin-echo  techniques,  but  this  requires 
that  laser  intensity  fluctuations  are  small.  Likewise,  by  ad¬ 
justing  the  laser  polarization  to  null  the  AC  Stark  shift  from 
a  single  beam,  we  mitigated  qubit  decoherence  due  to  laser 
intensity  fluctuations.  For  the  benchmarking  measurements 
described  here  we  did  not  actively  stabilize  (i.e.,  noise  eat)  the 
laser  beam  intensity. 

The  off-resonant  laser  beam  frequency  was  detuned  from 
the  cycling  transition  (|f)  -A  \P3/2,mJ  =  3/2))  by  A R  ~ 
—63.8  GHz.  This  gives  detunings  of  +15.6  GHz  and 
—26.1  GHz,  respectively,  from  the  |j~)  -A  P3/2 ,mj  =  1/2) 
and  14,)  -A  \P-}/2-,mJ  =  —  1/2)  transitions  (Fig.  4).  Laser  beam 
waists  were  wz  ~  1 10/rm  in  the  vertical  (z-direction)  and  wx  ~ 
1  mm  in  the  horizontal  direction.  Here  we  define  the  waist  as 
the  distance  from  the  center  of  the  beam  over  which  the  elec¬ 
tric  field  intensity  decreases  by  1/e2  (i.e.  I(z)  ~  ). 

With  the  small  2.4°  incident  angle  each  beam  makes  with  re¬ 
spect  to  the  plane  of  the  crystal,  this  provided  a  uniform  elec¬ 
tric  field  with  <  10%  intensity  variation  across  ion  crystal  ar¬ 
rays  with  N  <  250. 

The  ODF  laser  beams  were  linearly  polarized  at  nonzero 
angles  with  respect  to  the  7-axis.  Let 

Eu(r,t)  =  Zu Ev cos  [ku  -r- tout ) 

— .  ,  „  _  {  (7) 

EL(r,t)  =  £lElcos  lkL  ■  r  —  (i)Lt  J 

denote  the  electric  fields  of  the  upper  and  lower  ODF  beams. 
If  <])„  is  the  angle  of  the  laser  beam  electric-field  polarization 
with  respect  to  vertical  polarization  (e -x  =  0),  then  the  AC 
Stark  shift  of  the  qubit  states  when  illuminated  by  a  single 
beam  can  be  written 

At,acii  =Atcos2(c|)p)+Btsiir  (t>p)  ^ 

A|.,aCss  =  A^cos2  (c|)p)  +B|Sin2  (<f»p) 

where  Ay(A^)  is  the  Stark  shift  of  the  |t)(|i))  state  for  a  7t- 
polarized  beam  (e  parallel  to  the  z-axis)  and  B(B)  is  the  Stark 
shift  of  the  |t)(ID)  state  for  a  o-polarized  beam  (e  perpendic¬ 
ular  to  the  z-axis).  (Here  we  neglect  the  small  a  polarization 
(°=  sin 2. 4°)  that  exists  when  (f+  =  0.)  The  Stark  shift  of  the 
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qubit  transition  is 

=  (At  -  A;)  cos2  (<j)p)  +  (Bt  -  B±)  sin2  (<\>p)  .  (9) 


If  Af  A|  and  B- |-  —  B  have  opposite  signs,  there  is  an  angle 
which  makes  Aacss  =  0.  For  a  laser  detuning  of  Ar  =  —63.8 
GHz,  Aacss  =  0  at  4>p  ~  ±65°. 

With  Aacss  =  0  for  each  ODF  laser  beam,  we  exploit  the 
freedom  to  choose  their  polarization  in  order  to  obtain  a  state- 
dependent  force.  Specifically,  we  choose  Ejj  to  have  a  polar¬ 
ization  given  by  (f>„u  =  65°  and  /:/  to  have  a  polarization  given 
by  §pj  =  —65°.  In  this  case  the  interference  term  in  the  ex¬ 
pression  for  the  electric  field  intensity  ^Eu  +  E/  j  produces  a 

polarization  gradient  which  results  in  spatially  dependent  AC 
Stark  shifts 

(A-f  cos2  (typ)  - fit sin2  (<|)p))  2cos (8 k-z- Hr!)  n m 

(A;  cos2  (0P) -B;  sin2  (§p))2cos(Sk- z~  HR*)  ' 


for  the  qubit  levels.  Here  6k 


=  2k  sin 


is  the 


wave  vector  difference  between  the  two  ODF  laser  beams, 
Br  =  ®£/  —  is  the  ODF  beat  note,  and  (j)p  =  |<|)pjM  |  =  |(j)/,./|. 
The  spatially  dependent  AC  Stark  shift  produces  a  state- 
dependent  force  F-M(z,f)  =  F0^.i sin  (6k- z  —  /upt)  where 


Fot  =  -2  5 k  (At  cos2  (<|)p)  -  Bt  sin2  (<>p)) 

Foi  =  -2  5k  (A;  cos2  (<>p)  -B;sin2  (<>p)) . 

In  general  Fot  ±  ~FoV  We  operate  at  Ar  =  —63.8  GHz  where 
for  Aacss  =  0  we  also  obtain  Fa^  =  =  F0 

For  a  given  (\>pu ,  §pj ,  and  we  use  straightforward  atomic 
physics  along  with  well  known  values  for  the  energy  levels 
and  matrix  elements  of  9Be+  to  calculate  F0  as  a  function  of 
the  electric  field  intensity  Ir  =  |E/|2  =  \Eu\2  at  the  cen¬ 

ter  of  the  laser  beams.  For  0/„>  =  4.8°  and  Ir  =  1  W/cm2  , 
F„  =  1.4  x  10-23  N. 

Stronger  forces  can  be  generated  after  experimental  mod¬ 
ification  to  our  apparatus  to  permit  Qr  «  35°  .  At  this  an¬ 
gle  larger  detunings  ,ur  —  C0i  are  required  to  satisfy  Eq.  6  in 
Methods.  With  our  definition  Hi  =  /,  'Li<  jJi.j^]^Zjr  the  in¬ 
teraction  strength  between  two  ions  i  and  j  is  Jij/N.  Con¬ 
sider  for  example  N  =  217  ions,  cor  =  27t  x  45.6  kHz,  pow¬ 
ers  of  Ir  =  20  mW/beam  (12.5  W/cm2)  and  Qr  =  35°.  The 
spin-motion  entanglement  constraint  is  satisfied  by  a  detun¬ 
ing  of  /Jr  —  ©i  =  2tc  x  100  kHz.  In  this  case  we  obtain 
Jij/(2nN)  ~  (560Hz)  (do/djj)11 ,  where  do  ~  20^m  is  the 
typical  nearest  neighbor  separation. 


HI.  WAVEFRONT  ALIGNMENT 

The  ODF  laser  beams  produce  a  one  dimensional  (ID)  op¬ 
tical  lattice  characterized  by  the  effective  wave  vector  8 k  and 
beat  note  ,ur.  In  Sec.  II  we  assumed  that  6k  ||  z,  or  equivalently 
that  the  wavefronts  of  the  lattice  were  aligned  perpendicular 
to  the  z-axis  (magnetic  field  axis).  If  the  wavefronts  are  not 
normal  to  the  f-axis  as  sketched  in  Fig.  6,  then  the  time  de¬ 
pendence  of  the  ODF  seen  by  an  ion  in  the  rotating  frame 


Figure  6.  (a)  Sketch  of  the  ODF  laser  beam  geometry  used  to  gen¬ 
erate  the  one  dimensional  (ID)  traveling  optical  lattice,  (b)  Sketch 
of  the  ID  optical  lattice  wave  fronts  (red  lines).  These  wave  fronts 
need  to  be  aligned  with  with  the  ion  planar  array  (represented  by  the 
blue  dots).  Here  Xr  =  2tt/|AE|  ss  3.7/rm  and  Qerr  denotes  the  angle 
of  misalignment.  400  /rm  is  the  typical  array  diameter  for  N  ~  200 
ions.  With  the  wavefront  alignment  technique  discussed  in  the  text 
we  obtain  0err  <  0.05°. 


depends  on  the  (x,y)  position  of  the  ion.  This  complicates  the 
effective  spin-spin  interactions  generated  by  the  ODF  but  can 
be  adequately  mitigated  by  careful  alignment. 

Alignment  of  the  ODF  laser  beams  is  obtained  by  a  tech¬ 
nique  that  makes  top-view  images  (images  of  the  ion  reso¬ 
nance  fluorescence  scattered  along  the  magnetic  field)  from 
a  single  plane  of  ions  sensitive  to  ODF  wave  front  misalign¬ 
ment.  For  this  measurement  we  set  /jr  =  0  (stationary  ID  lat¬ 
tice)  and  detune  the  frequency  of  the  ODF  laser  beams  approx¬ 
imately  0.5  GHz  below  the  |f)  — >  \2P3/2  mJ  =  +3/2)  Doppler 
cooling  transition.  This  small  detuning  generates  sufficiently 


Figure  7.  Top-view  image  of  the  spatially  inhomogeneous  fluores¬ 
cence  from  a  single  ion  plane  produced  by  the  AC  Stark  shift  from  a 
static  (/ ir  =  0)  ODF  lattice  with  misaligned  wave  fronts.  Dark  bands 
are  regions  of  high  standing  wave  electric  field  intensity  (parallel 
to  the  dashed  yellow  line).  The  bright  horizontal  feature  bisecting 
the  center  of  the  image  is  fluorescence  from  the  weak  Doppler  laser 
cooling  beam  directed  perpendicular  to  the  magnetic  field.  The  im¬ 
age  was  obtained  by  subtracting  a  background  image  with  the  ODF 
beams  off  (Is  integration). 
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Figure  8.  (a)  Sequence  used  to  measure  spin  precession  due  to  the 
effective  mean  field  generated  by  the  engineered  Ising  interaction. 
The  spins  are  initially  prepared  in  |t).  The  first  pulse  rotates  the 
spins  by  an  angle  0]  about  an  axis  (defined  to  be  the  i-axis  in  the 
rotating  frame)  in  the  equatorial  plane.  The  subsequent  71  and  7t/2 
pulses  occur  about  the  y-axis.  (b)  Evolution  of  a  single  spin  prior  to 
the  application  of  the  spin-echo  7t-pulse. 


large  AC  Stark  shifts  on  the  cooling  transition  to  measurably 
change  the  ion  scatter  rate  from  the  Doppler  cooling  laser. 
With  the  Doppler  cooling  laser  turned  on  and  the  ODF  beams 
turned  off,  we  observe  a  spatially  uniform,  time-averaged  im¬ 
age  of  a  rotating  planar  crystal.  With  the  ODF  beams  on,  ions 
located  in  regions  of  high  electric  field  intensity  at  the  anti¬ 
nodes  of  the  optical  lattice  are  Stark-shifted  out  of  resonance 
with  the  Doppler  cooling  laser.  This  is  the  cause  of  the  dark 
bands  in  the  top-view  image  shown  in  Fig.  7.  We  adjust  the 
ODF  laser  beams  based  on  this  real-time  imaging  to  optimize 
their  alignment.  Improved  alignment  is  indicated  by  a  fringe 
pattern  of  longer  wavelength.  With  this  technique  we  have 
aligned  the  ODF  wave  fronts  with  the  planar  array  to  better 
than  ~  0.05°. 

The  image  in  Fig.  7  is  typical  of  what  we  obtain  with  1  s 
integration  duration.  This  indicates  the  ID  lattice  was  stable 
during  the  integration  period  and  shows  the  phase  stability  of 
our  ID  lattice  of  better  than  1  s. 

We  note  that  direct  fluorescence  imaging  of  the  ID  lat¬ 
tice,  for  example  by  tuning  the  ODF  laser  resonant  with  the 
Doppler  cooling  transition,  is  not  viable.  Even  at  low  powers, 
resonantly  scattered  photons  across  the  large  horizontal  waist 
of  the  ODF  beams  apply  a  large  torque,  causing  the  rotation 
frequency  and  radius  of  the  array  to  rapidly  change,  typically 
driving  the  ions  into  orbits  of  very  large  radius. 

We  have  also  used  phase  coherent  Doppler  velocimetry  to 
improve  the  ODF  wavefront  alignment  [45].  But  the  top-view 
imaging  technique  discussed  here  and  shown  in  Fig.  7  pro¬ 
vides  more  information  on  the  angle  and  direction  of  mis¬ 
alignment,  which  greatly  improves  the  ODF-crystal  alignment 
process. 


IV.  MODELING  MEAN  FIELD  SPIN  PRECESSION 


For  any  spin 


=  0,  which  implies  ( dj  )  =  con¬ 


stant  under  application  of  ///.  The  Hamiltonian  Ilf1  describ¬ 
ing  the  mean  field  response  of  spin  j  to  the  application  /?/  is 


obtained  by  expanding  Hi  to  first  order  in  8 o  ■  =  dj  — 


HfF  =  ^'£i<jJij((6zi}8c)  + 

LUn  £  -M6;  x8o;. 

U¥j 

Defining  Bj  =  ^  Ji,j  (d?),  HfF  can  be  written 

Hr=tsM/ 2. 

j=  1 


(12) 


(13) 


The  mean  field  Heisenberg  equations  of  motion  for  spin  Ga. 
is  obtained  from  the  commutator  of  d/;  and  Ilf1.  However, 

because  8 a  j  =  dy  —  (of)  and  (dj()  is  a  constant,  it  is  clear 
that  Eq.  13  describes  spin  precession  about  the  z-axis  at  fre¬ 
quency  Bk.  Our  observable  is  the  spin  precession  (in  excess 
of  ordinary  Larmor  precession)  averaged  over  all  the  spins 

1  x^N  z5 
uLk=\Bk- 

We  use  the  spin-echo  sequence  in  Fig.  8  to  measure  a  pre¬ 
cession  proportional  to  the  expectation  value  of  the  spin  pro¬ 
jection  along  z.  ((d?)).  The  spin-echo  sequence  minimizes 
contributions  to  spin  precession  that  are  not  oc(d?).  Specifi¬ 
cally,  the  spin  echo  cancels  a  constant  spin  precession  inde¬ 
pendent  of  (dj)  (e.g.,  due  to  slow  uncontrolled  magnetic  field 
fluctuations),  but  precession  (dy)  coherently  adds  in  the  two 
arms  of  the  sequence. 

The  first  pulse  sets  (dy)  =  cos  0i .  The  interaction  ///  is  then 
applied  by  turning  on  the  ODF  laser  beams  for  a  period  Xarm. 
During  this  interval  (dy)  is  a  constant  and  mean  field  theory 
predicts  an  average  spin  precession  angle  of  2/cos(0i)  •  Xarm. 
The  7t-pulse  changes  (dy)  — >  —  cos(0i)  and  the  spin  preces¬ 
sion  angle  2/cos(0i)  -Xarm  to  —  2/cos(0i)  • xarm ■  At  the  end 
of  the  second  Hi  interaction  of  duration  Xarm,  the  total  preces¬ 
sion  angle  of  the  spins  is  — 2/cos(0i)  -2xarm.  The  final  n/2- 
pulse  is  about  an  axis  shifted  by  90°  from  the  first  0i -pulse. 
This  pulse  converts  precession  out  of  the  initial  yz  plane  into 
excursions  above  or  below  the  equatorial  plane  of  the  Bloch 
sphere,  which  we  measure. 

The  evolution  operator  of  the  measurement  sequence  Useq 
is  obtained  in  a  straight  forward  manner  from  the  individual 
evolution  operators  from  each  segment  of  the  sequence. 


Useq  =k(Si\)-  U{HfF) -R(y,n)-U  (HfF)  -R(x,Q  i)  . 

(14) 

.  R(y,n)  = 


Here  7?(x,0i)  = 


cos(0i/2)  — fsin(0i/2) 

— / sin  (0i  /2)  cos  (0i /2) 


'0-1' 
1  0 

and  R  (y,  f )  =  4 

'  1  -1  ' 
1  1 

lution  is  given  by 

.  The  mean  field  evo- 


0  (HfF)  = 


exp  (  iJ  (dy)  xann)  0 

0  exp  (iJ (dj)  xarm) 


(15) 


where  (az)  =  cos(0i)  in  the  first  arm  of  the  spin-echo  se¬ 
quence  and  (d,)  =  —  cos(0i)  in  the  second  arm.  At  the  end 


10 


Figure  9.  We  determine  Ir  from  the  AC  Stark  shift  of  the  qubit  tran¬ 
sition  as  a  function  of  ODF  laser  beam  polarization  angle  §p.  The 
triangles  and  filled  circles  are  separate  measurements  done  on  the 
upper  and  lower  ODF  beams.  Ir  is  determined  independently  for 
each  ODF  beam  by  a  fit  to  Eq.  9.  The  error  bars  are  smaller  than  the 
plot  points. 

of  the  sequence  we  detect  the  |f)  state  probability.  Explicit 
computation  gives 

|(t  \Useq\  t)|2  =  \  {1  +  sin (0i)  -sin [2/cos (0i)  -2xarm]}  . 

(16) 

V.  OPTICAL  DIPOLE  FORCE  LASER  INTENSITY 
CALIBRATION 

The  spin-spin  couplings  (/)  measured  here  depend  on  the 
square  of  the  ODF  electric  field  intensity  Ir.  Therefore,  care¬ 
ful  calibration  of  the  electric  field  intensity  was  important  in 
benchmarking  the  strength  of  the  interactions  with  mean-field 
predictions.  We  separately  determined  Ir  for  each  beam  by 
measuring  the  AC  Stark  shift  of  the  qubit  transition  as  a  func¬ 
tion  of  the  laser  beam  polarization  angle  §p,  as  shown  in  Fig. 
9.  The  qubit  transition  frequency  was  measured  by  fitting  for 
the  center  frequency  of  a  Rabi-resonance  profile.  The  po¬ 
larization  angle  §p  was  varied  by  rotating  a  A/2  plate.  The 
Stark  shift  measurements  were  fit  to  Eq.  9,  which  provided 
values  for  A  and  B  (Stark  shifts  for  7t-polarization  and  o- 
polarization).  Values  for  A  and  B  were  then  used  with  atomic 
physics  calculations  to  determine  Ir,  the  electric  field  intensity 
of  the  ODF  laser  beams,  with  a  fractional  uncertainty  of  ~  5%. 
Frequent  intensity  calibration  measurements  were  taken  dur¬ 
ing  a  benchmarking  run.  Slow  drifts  to  the  laser  intensity  be¬ 
tween  calibrations  added  another  5%  uncertainty,  which  we 
add  in  quadrature  to  the  fitted  uncertainty. 

VI.  SPONTANEOUS  EMISSION 

Decoherence  due  to  spontaneous  emission  has  been  well 
studied  in  this  system  [41].  The  qubit  levels  are  closed  under 
spontaneous  light  scattering;  that  is,  spontaneous  light  scat¬ 
tering  does  not  optically  pump  the  ion  to  a  different  ground 


state  level  outside  of  the  two  qubit  levels.  We  measure 
spin  precession  to  benchmark  the  Ising  interaction  couplings. 
Equivalently  we  measure  the  evolution  of  off-diagonal  coher¬ 
ences  between  the  |f)  and  |)  levels.  Spontaneous-emission- 
induced  decay  of  these  coherences  is  accurately  modeled  by 
Eq.  (8)  of  Ref.  [41],  We  add  this  time  dependence  to  the  uni¬ 
tary  evolution  U  ( H ^F)  discussed  in  Sec.  IV.  This  results  in 
the  expression 

p(|t))  =  \  (1  +  exp (— r- 2xarm) sin (0i) sin (2/ cos (0i)  •  2xorm))  , 

(17) 

where  P(|t))  is  the  probability  of  detecting  a  spin  in  ||) 
at  the  end  of  the  spin-echo  sequence.  The  only  difference 
with  Eq.  16  is  the  factor  exp(— T  •  2 Tam)  due  to  the  decay  of 
the  off-diagonal  elements  of  the  density  matrix.  Here  T  ac¬ 
counts  for  decoherence  due  to  spontaneous  emission.  From 
[41],  r  =  i  (rRam+rel)  has  contributions  from  both  Raman 
scattering  and  elastic  Rayleigh  scattering.  A  straightforward 
atomic  physics  calculation  relates  T  to  the  ODF  laser  beam 
intensity  Ir.  In  fits  of  Eq.  17  to  the  spin-precession  measure¬ 
ments,  we  fix  r  at  the  value  determined  by  the  laser  intensity 
calibrations.  For  A r  =  —63.8  GHz,  (j)p  =  ±65.3°,  and  Ir=  1 
W/cm2,  we  calculate  T  =  82  s~  1 . 

With  our  present  set-up  we  can  achieve  /^>  T  for  detunings 
| /jr  —  ©i  |  <  10  kHz.  For  these  detunings  it  will  be  possible  to 
simulate  quantum  effects  beyond  mean  field  theory  such  as 
spin  squeezing  and  spin-depolarization  due  to  many-body  in¬ 
teractions.  For  example  we  calculate  that  with  4  mW/beam 
and  an  ODF  beatnote  detuning  of  \/jr  —  (Oi  |  =2  kHz  ,  J  ~ 

271  x  0.5  kHz  and  T /J  ~  0.06.  The  potential  spin  squeezing 
due  to  this  interaction  is  5  dB,  limited  by  spontaneous  emis¬ 
sion.  For  many  simulations  it  will  be  desirable  to  achieve 
J/N  3>  r  where  /  is  a  nearest-neighbor  coupling  strength. 

(J  ~  /  for  small  detunings  \/jr  —  <Mi  | .)  The  most  straight  for¬ 
ward  strategy  to  achieve  this  condition  in  our  set-up  appears 
to  be  to  increase  the  angle  Or  between  the  Raman  beams.  The 
ratio  J/r  scales  as  sin2  (Qr/2).  Therefore  decoherence  due  to 
spontaneous  emission  can  be  dramatically  reduced  with  an  or¬ 
der  of  magnitude  increase  in  Qr.  We  note  that  a  large  increase 
in  Qr  will  require  implementation  of  a  sub-Doppler  cooling 
scheme  to  remain  in  the  Lamb-Dicke  regime  (see  Sec.  VII). 

Different  ODF  laser  detunings  can  likely  help  reduce  the 
impact  of  spontaneous  emission  (with  the  complication  that 
Ff  /  Fj,  but  very  large  laser  detunings  A  obtained  by  tun¬ 
ing  the  ODF  laser  beam  frequencies  outside  the  P\/2  —P3/2 
manifold  appear  unlikely  to  help.  This  is  because  both  the 
interaction  strength  and  spontaneous  emission  scale  as  1/A2 
for  our  qubit,  which  is  not  a  clock  transition.  For  trapped-ion 
experiments  in  low  magnetic  field,  the  impact  of  spontaneous 
emission  is  typically  minimized  by  tuning  the  ODF  laser  beam 
frequencies  between  the  P\/ 2  and  P3/2  manifolds.  Because  the 
9Be+  P-state  fine  structure  (197  GHz)  is  comparable  to  the 
Zeeman  splittings  for  our  magnetic  field  (4.5  T),  this  approach 
does  not  significantly  help.  However,  it  can  help  with  heavier 
mass  ions  (e.g.,  Mg+).  Reference  [37]  discusses  a  number  of 
strategies  that  can  minimize  the  impact  of  spontaneous  emis¬ 
sion  when  the  P-state  fine  structure  is  large  compared  to  the 
Zeeman  interaction  energy. 
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VII.  LAMB-DICKE  CONFINEMENT 

Throughout  this  manuscript  we  implicitly  assume 
that  the  applied  state-dependent  force,  = 

cos  (8£  •  z  —  /jRt),  is  constant  across  the  spatial  ex¬ 
tent  of  an  ion’s  wave  function.  Specifically,  for  a  planar  array 
located  at  z  =  0,  we  assume  the  force  on  any  ion  is  given  by 

Ft,f  (z  =  0,f)  =  F0 ucos  (/uRt)  .  (18) 

The  extent  to  which  this  is  true  is  quantified  by  the  parameter 

h =  8/\  ■  Zrm\.i  (19) 

where  zrms,i  is  the  root  mean  square  (rms)  axial  extent  of  the 
wave  function  of  ion  i  and  8k  =  2tt/A ,r,  where  Xr  =  3.7 /.mi  for 
0r  =  4.8°.  We  note  that  r\ind,i  is  a  Lamb-Dicke  confinement 
parameter  for  a  individual  ion  (ind),  not  to  be  confused  with 
the  usual  Lamb-Dicke  parameter  that  is  defined  in  terms  of 
the  ground-state  wave  function.  In  the  limit  that  the  Coulomb 
interaction  energy  between  ions  is  a  small  perturbation  to  the 
axial  potential  of  the  trap,  we  can  think  of  each  ion  as  a  single 
ion  confined  in  the  external  trap  potential,  and  approximate 

ZrmsJ  —  '  V^n+l,  where  M  is  the  mass  of  an  indi¬ 

vidual  ion.  Recently  we  have  completed  careful  axial  tem¬ 
perature  measurements  for  single-plane  ion  crystals  [30].  We 
find  ~  1  mK  for  the  axial  COM  mode  and  ~  0.4  mK  for 
the  higher  order  transverse  modes.  For  (Dz  =  2n  ■  795  kHz 
and  conservatively  assuming  T  =  1  mK,  we  calculate  n  ~  26, 
Zrms,i  -  190  nm,  and  V\mdl  ~  0.32. 

For  the  benchmarking  measurements  described  here  we 
typically  set  the  rotation  frequency  of  the  ion  array  about 
~  0.5  kHz  below  the  rotation  frequency  of  the  1  f>2  plane 
transition,  determined  experimentally  from  side-view  images 
[21],  In  this  case,  the  spectrum  of  the  transverse  (axial)  modes 
is  broad,  and  we  underestimate  Zrms,i  in  the  above  analysis.  An 
improved  estimate  of  Zrms,i  is  obtained  by  summing  the  con¬ 
tributions  from  all  of  the  transverse  modes  m 

(  h  \1//2 

Zrmsj  =  (  Y.  ( bijn  f  ZZTJ (2 nm  +  1)  )  ,  (20) 

\m  2M®m  ) 

where  hm  ~  kpT /Tmln  is  the  mean  thermal  occupation  of 
mode  m.  We  assume  here  that  every  mode  is  characterized  by 
the  same  temperature  T.  For  A  =  217  and  co.  =  2n  x  795  kHz, 
the  rotation  frequency  of  the  1  f>2  plane  transition  is  0),  = 
2n  x  46.1  kHz.  For  (Or  =  2n  x  45.6  kHz,  the  spectrum  of  the 
axial  modes  ranges  from  795  kHz  down  to  224  kHz.  With 
Eq.  20  we  calculate  zrrm.i  —  520  nm  in  the  center  of  the  array, 
decreasing  to  Zrms,i  —  250  nm  at  the  array  edge,  corresponding 
to  q i„dj  ~  0.89  for  i  near  the  array  center  and  v\ind.,  =  0.42  for 
i  near  the  array  boundary. 

The  wavefront  alignment  discussed  in  a  previous  section 
can  be  viewed  as  achieving  a  type  of  Lamb-Dicke  confine¬ 
ment.  Let  0 denote  the  angle  of  misalignment  between  the 
planar  array  and  the  ODF  ID  lattice  (Fig.  6).  Relative  to  the 
ODF  lattice  wavefronts,  the  rotation  of  the  array  produces  a 
time-dependent  shift  in  the  axial  position  of  an  ion  that  can 


be  written  as  z(t)  =  Rsin(0f,T)sin(GV+(p).  Here  R  is  the 
distance  of  the  ion  from  the  center  of  the  array,  and  cp  is  deter¬ 
mined  by  the  azimuthal  position  of  the  ion  in  the  array.  The 
force  on  this  ion  is  then 

#r,4  (0  =  F°trl cos  (8fe/?sin0err  •  sin(GV  +9)  ~FRf)  •  (21) 

ForEq.  21  to  approximate  Eq.  18,  we  desire  8k-2R/>sin0m.  < 
1,  where  Rp  is  the  array  radius.  For  Rp  ~  200  /./m  (typical  for 
N  =  200)  and  0,.,,-  ~  0.05°  (see  Wavefront  alignment  section), 
we  calculate  8k-2RpsinQerr  ss  0.6. 


VIII.  ION  LATTICE  CONFIGURATION  AT 

EQUILIBRIUM  AND  TRANSVERSE  NORMAL  MODES 

The  ion  equilibrium  positions  r,-  at  zero  temperature  are  cal¬ 
culated  by  minimizing  the  Euler-Lagrange  action  for  the  Pen¬ 
ning  trap  potentials  and  the  constraint  that  the  ions  lie  in  a 
plane  (at  z  =  0).  The  solution  is  a  triangular  lattice  with  a 
lattice  constant  that  increases  as  one  moves  radially  outward, 
that  has  a  smooth  unfaceted  edge  and  exhibits  a  degradation  in 
orientational  order  near  the  crystal  perimeter.  The  transverse 
(along  z)  phonon  modes  (co,„,  bm)  are  obtained  by  Taylor  ex¬ 
pansion  of  the  potential  about  the  ion  equilibrium  positions  r,. 
In-plane  modes  (along  x,y )  can  also  be  calculated,  by  solv¬ 
ing  a  quadratic  eigenvalue  problem  (due  to  inclusion  of  the 
centrifugal  and  Coriolis  forces).  The  details  of  the  transverse 
mode  calculation  is  discussed  in  this  section.  The  problem  has 
also  been  solved  for  longitudinal  and  transverse  modes  in  ID 
(e.g.,  [46,  47]). 

In  general,  the  Lagrangian  for  a  collection  of  N  ions  with 
charge  q  and  mass  M  in  an  electromagnetic  field  6,  —  A(r,)  is 

L  =  T-V=  f  {  l-MXf  -  q  (<>,■  | ,  (22) 

where  f,-  =  (r,.  0,,z,j  is  the  coordinate  of  ion  i.  In  a  Penning 
trap,  the  field  consists  of  a  uniform  magnetic  field  in  the 
f-dircctioii  (B  =  Bqz)  and  a  harmonic  trapping  (anti-trapping) 
electric  potential  in  the  ’-direction  (f-direction)  with 
frequency  coi .  Lastly,  an  additional  external  time-dependent 
electric  quadmpole  potential  (the  “rotating  wall”,  amplitude 
Vwan  at  the  ions)  is  applied  to  control  the  ion’s  rotation 
frequency  in  the  trap  [48],  Thus,  the  scalar  potential  for  ion  i 
is 

ytyirap.i  T  qtywalt.i  ytycoulombj 
=  (z2i  -  rj /2)  +  Vwaiirf  cos 2(0,-  +  C 0rt )  + 

where  the  ion-ion  separation  is  fjd  =  —  r and  k—  1  /47t£o- 

Since  B  is  uniform,  the  vector  potential  energy  is 

qA(?i)  =  -|q-  X  B .  (24) 

The  position  of  the  ions  at  equilibrium  is  calculated  as  fol¬ 
lows.  We  move  to  a  rotating  frame  where  the  ions’  coordi¬ 
nates  are  stationary  by  using  the  coordinate  transformation 
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r'j  =  (r-,0'-,z-)  =  (7,,  0  +  GVf,z,-),  a  counterclockwise  rotation. 
The  ions’  equilibrium  positions  r-  can  be  found  by  solving  the 
transformed  Euler-Lagrange  equations  L' .  To  seed  the  numer¬ 
ical  solution,  we  supply  an  initial  guess  for  the  2D  crystal:  a 
regular,  triangular  lattice.  Consistent  convergence  requires  ion 
numbers  N  corresponding  to  closed  shells  (e.g.,  N  =  127  has 
six  closed  shells).  We  find  ion  equilibrium  positions  that  devi¬ 
ate  from  a  perfect  triangular  lattice  near  the  crystal  periphery 
and  which  have  an  overall  ellipticity  due  to  Vwau. 

Given  7/  the  crystal’s  transverse  eigenmodes  can  be  calcu¬ 
lated  by  Taylor  expansion  of  the  potential  about  the  equilib¬ 
rium  positions  7/  The  Lagrangian  is 

C  =  1£at^'-I  £  Wyz'z'  (25) 

!=  1  Z  1,7=1 

where  zj  is  the  axial  displacement  of  ;-th  ion  and  K\j  =  Kjj  is 
the  symmetric  stiffness  matrix  evaluated  for  the  equilibrium 
configuration  as 


Kij  = 


kq 2 


kq2 


-1  \^n 


i  =  j,n  ^  i 

i^j 


(26) 


By  minimizing  the  action  8  j  dtL'  =  0  with  respect  to  the  ro¬ 
tating  frame  axial  coordinates  zj,  we  obtain  N  equations  of 
motion 


?i+t^j  =  0,fori=l,2,...,N.  (27) 

Pi  M 

Following  standard  normal  mode  analysis,  the  solution  is  ob¬ 
tained  by  calculating  the  eigenmodes  of  the  matrix  Kjj/M. 
The  result  is  N  eigenvalues  to2  ,  for  m  =  1,2,  ...,1V  (with  cor¬ 
responding  frequencies  co). 


IX.  LIMITS  TO  THE  VALIDITY  OF  MEAN  FIELD 
THEORY 

Here  we  calculate  the  domain  over  which  a  mean  field  the¬ 
ory  (MF)  treatment  of  Jr  is  valid.  In  the  special  case  of  uni¬ 
form  Ising  coupling  (e.g.,  0  <  /jr  —  CDj  <C  coi  —  CO2),  we  have 

f  (g^)  (la)/2)  =  08) 

where  we  have  neglected  a  constant  offset  and  /,  =  }_'v_ ,  d-/2 
is  the  z-component  of  the  total  composite  spin  of  the  system. 
The  interaction  strength  %  is  independent  of  ion-ion  separation 
and  is  given  by 


ll _ 1__ 

KIM  n ^  —  co2 


(29) 


We  start  each  experiment  with  N  spins  all  prepared  in  the 
|t)  state.  Here  we  use  the  composite  spin  picture  where  this 
state  is  labeled  \J  =  N/2,  Mj  =  N/2).  (Formally  \J,Mj )  la¬ 
bels  a  state  that  is  an  eigenstate  of  J2  and  J7  with  eigenvalues 

7(7+  1)  and  My,  where  7  =  £,-d;/ 2  is  the  total  system  spin.) 
The  first  pulse  of  the  spin-precession  measurement  sequence 
rotates  the  composite  Bloch  vector  by  an  angle  0/.  After  the 
rotation  by  0i  the  state  can  be  written  as 

N/2 

l¥>=  E  C(J,Mj)\J,Mj),  (30) 

Mj=—N/2 


where  the  coefficients  C(7,  Mj)  are  significantly  non-zero  for 
a  small  range  AM/  centered  on  m\°\  where  M «  ^cosO/. 
For  a  coherent  spin  state,  AM/  <  \/N . 

We  want  to  establish  that  a  7r  interaction  looks  like 

precession  °e  m\°\  at  least  over  short  durations.  A  precession 
by  an  angle  <f>  about  the  z-axis  is  obtained  with  the  operation 

g-'Wz 

e~ U  |¥>  =  EC(^  Mj)e ~i¥z \J,  Mj)  ■  (31) 

Mj 

We  measure  expectation  values  of  the  total  system  spin  7  .  It 
is  sufficient  therefore  to  consider  how  these  expectation 
values  transform  under  the  rotation  el0j?, 

(tj/|7v  |\|/)  ->  (\|/|  e^fzJye~^Jz  |\|/>  ■  (32) 

For  simplicity,  only  matrix  elements  of  the  total  spin  in  the 
y-direction  Jy  are  considered.  Identical  expressions  are  also 
obtained  for  Jx. 

(¥lA-l¥) 

=  c(j >  m'j)*C(J,  Mj)  ( J >  Mj  \fy \J,  MJ ) 

=  E  Mj  {C(7,  Mj  -  1)*C(7,  Mj)  (7,  Mj-  1 1  Jy  |7,  Mj) 
+C(7,  Mj  +  1  )*C(7,  Mj)  (7,  Mj  +  1 1  Jy  |7,  Mj)  )  . 

(33) 

The  double  sum  is  eliminated  by  making  use  of  the  property 
that  Jy  has  only  non-zero  matrix  elements  between  states 
with  Mj  —  Mj  =  ±1.  C(7,  Mj  ±  1)  can  be  defined  to  be  0  if 
Mj  —  1  =  — N/2  —  1  or  Mj  +  1  =  N/2  +  1.  Similarly 

(¥|  eli>J'zJye~l^z  |\|/) 

=  c(7,  M'jfCiJ,  Mj)eiM>e-iM *  ( 7 ,  m'j  |7y|  7,  Mj) 

=  E mj  {C(7,  Mj  -  1  )*C(7 ,  Mj)e~*  (7,  Mj  -  1 1  Jy  |7,  Mj) 
+C(J,  Mj  +  1)*C(7,  Mj)e*  (7,  Mj  +  l\  fy  |7,  Mj) }  . 

(34) 

This  is  to  be  compared  with  the  same  matrix  element  under 
the  Hi  =  j/fj'2  interaction. 


,  2 


=  Em'  Mj  c(7,  MjfC{J,  Mj)e &  K)  (7,  Mj \jy  |  7,  Mj 

{c(7,  My  -  1)*C(7,  Mj)e~*J*J  (-Mj  —  l)t  My  -  1  |7V  |7,  My) 

+C(7,  My  +  1  )*C(7,  My)e'T  (7,  My  +  1 1  7V  |7,  My)  } 


-My, My 

=  Em,, 


(35) 


X 
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Equation  35  approximates  a  rotation  about  the  z-axis  (that  is, 
approximates  Eq.  34)  if 

(2 Mj  -  1)  «  (2 Mj  +  1)  «  ■  2  Mf]  for  all  Mj  for 

which  C(J.  Mj)  is  significantly  non-zero.  This  is  satisfied  for 
short  periods  t  satisfying  jjjt  ■  2  (AM /  )  <C  1 .  For  the  initial 
coherent  spin  state,  A Mj  <  \/7V,  which  puts  a  limit 


(36) 


on  the  period  t  for  which  the  I:  interaction  can  be 
approximated  as  a  precession  ;V/(n  '.  The  precession 
frequency  predicted  by  the  /r  analysis, 

■  2 =  jjf  ■  cos 0  =  2% cos 0,  is  in  agreement  with  the 
precession  frequency  predicted  by  the  MF  analysis. 


Equation  36  provides  a  limit  on  the  period  for  which  the  MF 
analysis  is  valid.  We  confirm  this  limit  through  explicit  cal¬ 
culation  of  the  spin  precession  benchmarking  sequence  for  the 
Hi  =  jjjJ:  interaction  with  N  =  5,  50  and  100  spins.  The  exact 
calculation  is  possible  because  the  system  remains  in  the  sym¬ 
metric  subspace  (dimension  of  Hilbert  space  is  2N  ■  }  1).  Fig¬ 
ure  10(a)  shows  the  results  of  the  calculation  for  five  spins. 
Reasonable  agreement  between  the  MF  precession  formula 
(Eq.  17  with  T  =  0)  and  the  exact  calculation  is  obtained  for 
yt  =  0.2 ,  but  not  for  yt  =  0.8  and  1 .6,  as  expected  from  Eq.  36. 
Figures  10(b)  and  10(c)  compare  the  exact  calculation  with 
the  MF  formula  for  N  =  50  and  100  spins.  Excellent  agree¬ 
ment  between  the  MF  formula  and  the  exact  calculation  is  ob¬ 
tained  for  yt  =  0.2  and  0.8.  Some  differences  are  observed  at 
yj  =  1.6.  These  differences  decrease  as  N  increases. 


5  Particles 


50  Particles 


100  Particles 


Figure  10.  (7Z)  vs  initial  polar  angle  0]  for  an  initial  coherent  spin 
state  of  N  spins,  t  =  2zarm  in  the  legend  is  the  total  interaction  period. 
The  solid  line  is  the  exact  calculation.  The  open  circles  are  from  the 
MF  precession  formula  (Eq.  17,  with  T  =  0).  (a)  N  =  5.  (b)  N  =  50. 
(c)  N  =  100. 
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